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We review recent progress in the theory of neutron matter with particular emphasis 
on its superfluid properties. Results of quantum Monte Carlo calculations of simple 
and realistic models of uniform superfluid neutron gas are discussed along with those 
of neutrons interacting in a potential well chosen to approximate neutron-rich oxygen 
isotopes. The properties of dilute superfluid Fermi gases that may be produced in atom 
traps, and their relations with neutron matter, are illustrated. The density dependence 
of the effective interaction between neutrons, used to describe neutron-rich systems with 
the mean field approximations, is also discussed. 



1. INTRODUCTION 

Neutron matter is believed to be an unusual superfluid gas with positive pressure at 
all densities. It has been studied in the context of neutron stars for many years. Recent 
reviews can be found in references [1-3]. However, the subject is rapidly progressing; the 
large number of authors in this paper represents the breadth of recent interest. 

The interesting properties of neutron matter follow from the fact that the nn scattering 
length (a = —18 fm) has a large magnitude, so that the dimensionless parameter \akp\ 3> 1 
at densities > 10 _4 po ; where kp is the Fermi momentum given by kp = 37r 2 p. We use 
Po = 0.16 fm -3 to denote the density of nuclear matter in large nuclei. On the other 
hand, the effective range of the nn interaction (r nn ~ 2.8 fm) is not much smaller than 
the interparticle distance r = (3P/47T) 1 / 3 , in neutron matter at densities of interest. Hence 
the nn interaction can not be easily approximated by a zero range 5-function interaction. 
Realistic models of the nn interaction have a repulsive core and strong tensor component; 
thus simple perturbation theory expansions can not be easily used to describe neutron 
matter. In addition, due to the strong nn attraction in 1 Sq channel, the low-density 
neutron matter is believed to be an S— wave superfluid, although at higher densities it 
could be a P— wave superfluid [3]. 

Recently it has become possible to study cold dilute Fermi gases in atom traps [4]. In 
these gases the interatomic interaction can be approximated by a zero range potential 
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characterized by the scattering length a. The experiments have the ability to tune values 
of akp by using Feshbach resonances. These studies essentially probe superfluid gases like 
neutron matter at very low densities (r ^> 2.8 fm). 

Density-dependent effective interactions are commonly used in nuclear physics to de- 
velop mean-field approximations based on energy-density functionals. In absence of exper- 
imental data, results of many-body calculations of the equation of state of pure neutron 
matter [5,6] are often used to determine the effective interactions [7-9]. 

In section 2 we briefly describe the progress in quantum Monte Carlo calculations of 
the properties of dilute Fermi gases and neutron matter. Section 3 is devoted to quantum 
Monte Carlo calculations of neutrons in a potential well aimed to approximate neutron- 
rich oxygen isotopes. We discuss effective interactions in neutron matter in section 4, and 
conclude in section 5. 

2. UNIFORM SUPERFLUID GASES AND NEUTRON MATTER 

Recently it has become possible to calculate the energies and gaps of simple super- 
fluid Fermi gases ab initio from the interparticle interaction, using quantum Monte Carlo 
methods [10]. The total energy is calculated for N particles in a periodic box of volume 
N/p. The results for dilute (effective range <C r ) Fermi gas with akp = oo are shown in 
Fig.1. The unit of energy: 

e fg = ^--4 , (1) 

10 m 

corresponds to the energy per partiele of noninteracting Fermi gas. Dimensional argu- 
ments show that the energy and gap of dilute superfluid gases scale with EpQ when 
a — > oo. The pairing gap A is obtained from the difference between results for odd and 
even N (Fig.1). 

These studies of dilute superfluid Fermi gases have been extended to various values of 
akp in the domain — oo < (akp)' 1 < oo. Cubic periodic boxes have shell structure with 
shell closures at 2, 14, 38, ... for spin-| fermions. In dilute gases the pairing removes shell 
effects when the interaction is attractive and strong enough so that (akp)' 1 > —0.5. The 
gap energy defined as: 

A(N) =E(N) -^-[E(N -1) + E(N + 1)} , odd N , (2) 

does not depend significantly on N for (akp)' 1 > —1. At negative values of akp, the 
calculated values of A are in between the estimates obtained from BCS and Gorkov 
equations [10]. However, at large positive values of (a/c^) -1 , the gap energy is below the 
Gorkov estimate and approaches \Ei\/2, where Ei is the energy of the two-body ground 
state. Recall that when the interaction is attractive and a > 0, there must be one or more 
two-body bound states. 

As (akp)' 1 varies from — oo to +oo the character of dilute Fermi gases changes from BCS 
supernuids with weak pairing to diatomic Bose gases with Bose-Einstein condensation. 
This transition has been discussed in the literature [11,12]. 

The pairing in nuclei and neutron matter is much weaker than in the dilute Fermi 
gases with (akp)' 1 > —0.3. However, the energy of neutron matter appears to be close 
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Figure 1. GFMC results for the E(N) of dilute superfluid gas with ctfci? — > oo, in a 
periodic box. 

to the energy (~ 0A4Efg) predicted by quantum Monte Carlo calculations (Fig.1) for 
large over a wide density range. Most calculations of neutron matter with realistic 

interactions, carried out since 1970 [1], give E(p) ~ 0.5Efg{p) m the density range 0.01 
to 0.5po- The pairing gaps in neutron matter at p > O.lpo are sensitive to the details of 
the nn interaction. For example, the gaps predicted from only the scattering length and 
effective range are much larger than those with realistic models with repulsive cores [3]. 

Quantum Monte Carlo calculations of neutron matter with realistic interactions are 
much more difflcult than those for dilute gases. In order to treat the effects of the spin, 
tensor and spin-orbit components in the nn interaction exactly, one has to sum over all 
the 2 N spin states of the N interacting neutrons [13]. Such Green's function Monte Carlo 
(GFMC) calculations are limited to systems with no more than 16 neutrons at present. 
Methods to sample the spin space using auxiliary fields are also being developed. These 
auxiliary field diffusion Monte Carlo (AFDMC) calculations [14] are fairly accurate for 
pure neutron systems, and can be carried out for much larger values of N. 

In computations using periodic boxes the interparticle interaction beyond L/2, where 
L is the box length, has to be truncated or approximated. GFMC calculations [13] 
use truncated interactions and estimate the contribution of v (r > L/2) by approximate 
methods. It is small for p < po/2, but increases rapidly with density. The studies of pairing 
gaps in neutron matter have therefore focused on p < po/4. The present, preliminary, 
results are listed in Tab le I. 

The effect of shell closure at N = 14 is obvious in the results shown in Table I. The 
energy of the seventh pair, £(14) — £'(12) ~ 3.2 MeV, is much smaller than that of the 
eighth pair, £'(16) — £(14) ~ 26.8 MeV. Calculations with smaller statistical (Monte 
Carlo) error are necessary to study the dependence of A(N) on N, in addition to studies 
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Table 1 

Preliminary results of quantum Monte Carlo calculations of the energies of N neutrons at 
density po/4, in a cubic periodic box using the Argonne v' 8 interaction truncated at L/2, 
in MeV. 
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of A(N > 15) with the AFDMC method. Nevertheless the A(p /4) appears to be of order 
2.8(5) MeV. This value is in excellent agreement with that obtained using realistic bare 
nn interactions and the BCS equation, 2.6(1) MeV. This agreement is surprising since 
most estimates [3] suggested that medium polarization corrections to the bare interaction 
would suppress the gap by a factor of two or more even at a density of po/4- 

The BCS equation predicts a maximum gap of 3 MeV at 0.1p , and the disappearance 
of the 1 So pairing at ~ 0.7po [3]. It would be most interesting to verify these predictions 
with essentially exact ab initio calculations. 

3. Neutron Drops : Models of Neutron-Rich Nuclei 

Neutron drops bound in external wells were studied previously with various objectives 
[15-17]. We are presently exploring the possibility to model neutron-rich oxygen isotopes 
by neutron drops in a well meant to represent the effect of interactions with the protons. 
The model Hamiltonian for the isotope with N neutrons has the form: 



H = E + + E % + E v*k , (3) 

i=l,N \ Zm / i<j<N i<j<k<N 

where Vij and are realistic bare two-, and three-neutron interactions. The Argonne 
vis and Illinois-2 interactions are used in the present calculations and the single-neutron 
potential V(r) is chosen so that the density distribution of the eight-neutron ground state 
is close to that of protons in 16 O. It is given by: 

V ' W = l+expnf-3)/!.!] ' (4) 

in MeV and fm units. This well has only two bound states; a Os and a Op at ~ -12.1 and 
-4.8 MeV respectively. It can thus bind only 8 noninteracting neutrons. However, bound 
states of many more interacting neutrons exist in this well due to the nn attraction. 

The calculated energies of the low-energy bound states of up to 14 neutrons are com- 
pared with the experimental energies of oxygen isotopes in Fig.2. A constant of -19 MeV 
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Figure 2. The energies of neutron drops in the chosen external well are compared with 
the experimental energies of oxygen isotopes with the same number of neutrons. 



has been added to the calculated energies to match the energy of the 8-neutron drop 
with that of 16 0. The results indicate that neutron-rich (N > 8) oxygen isotopes can be 
well approximated by neutron drops in a Woods-Saxon well. However, the approxima- 
tion becomes less accurate for neutron-poor (N < 8) isotopes. Presumably the proton 
distribution in oxygen isotopes does not change much with N when N > 8, but when N 
becomes smaller than 8, the protons become less bound and their distribution changes 
significantly. 

In the traditional shell model, neutron-rich oxygen isotopes are approximated by N — 8 
neutrons interacting via an effective interaction in the sd-shell. In contrast the present 
method uses bare interactions among all the neutrons. The single-partiele potential V(r) 
is meant to represent only the np interaction effects. It can obviously be improved by 
adding a spin-orbit term to the V(r). 

In O the interactions between neutrons give only ~ 2/3 of the observed ^3/2—^5/2 spin- 
orbit splitting, and the three-neutron interaction contributes approximately half of the 
calculated splitting. The remaining third must come from the interactions with protons. 
The variation of the pairing gap with iV is being studied in more detail. 



4. The Effective Interaction in Neutron Matter 

Density-dependent effective interactions are often used to describe nuclei using the 
mean-field approximation. They are mostly used in the Hartree, Hartree-Fock or Hartree- 
Fock-Bogolubov equations. In principle, parts of such effective interactions which depend 
only on the density of matter, can be chosen such as to reproduce the E(p) and the pairing 
gap A(p) of uniform neutron matter. 
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The simplest Skyrme effective interaction between neutrons is of the form: 



v eff (R,r)=b(p(R)) 5(v). 



(5) 



Here R and r are the center of mass and relative separation between the two interacting 
neutrons. Commonly used Skyrme interactions have momentum dependent and spin- 
orbit terms which are omitted here for brevity. Their dependence on the proton density 
can also be neglected in the present discussion. In the Hartree-Fock approximation, the 
energy of neutron matter per neutron given by this v e ff is: 



The strength of the Skyrme interaction b(p) can thus be obtained from the E(p) of neutron 
matter. 

Most realistic calculations of neutron matter give E(p) ~ 0.5Efg(p) m the density 
range p = 0.01 to 0.1 fm~ 3 . This result is independent of the interaction model and 
the calculation details [1]. It simply follows from the fact that in this density range the 
|ö nn ^F| 1- This implies that: 



since E FG is proportional to p 2 / 3 . This low-density divergence of the effective nn interac- 
tion is also a consequence of the large magnitude of the nn scattering length. 

In some models [7] the neutron-proton effective interaction has a more severe (1/p) 
divergence due to the existence of the np bound state, the deuteron. However, this 
divergence is mostly ignored because matter containing neutrons and protons is unstable 
at low densities; it clusters into nuclei. In contrast the p -1 / 3 divergence of the effective 
nn interaction could have consequences on the structure of neutron-rich nuclei with a 
neutron halo and in the neutron gas in a neutron-star crust [1]. 

It is known that pairing is very important in determining the structure of neutron-rich 
nuclei. For example, nuclei with an even number of neutrons dominate the neutron- 
drip line. It is not possible to use delta-function effective interactions of finite strength 
to calculate pairing gaps without restricting the Hilbert space. In unrestricted space 
divergences occur. In stable nuclei reasonable values of the gaps may be obtained using 
the Skyrme effective interaction and truncating the Hilbert space to a single open shell 
[16]. However, in nuclei near the neutron drip line, shell structure becomes less dominant, 
and the problem becomes more difficult. Alternative approaches of using a local density 
approximation for the pairing fields [18] are also being pursued. 

5. Conclusions 

Advances in quantum Monte Carlo Methods in the new millennium have made possible 
essentially exact, ab initio calculations of the equations of state of low-density neutron 
matter. We are also developing methods to calculate the pairing gaps in low-density 
neutron matter, and in dilute Fermi gases from the bare interactions. All realistic models 
of nuclear forces predict similar properties of low-density neutron matter, thus these ab 



E(p) = E FG (p) + -b(p)p . 



(6) 



b(p) ~ — Epc(p) ~ ap 
p 



(7) 
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initio predictions are expected to have little model dependence. Ways to define effective 
interactions in neutron matter from the results of these ab initio calculations are being 
refined. They can presumably be used to study the structure of neutron-rich nuclei. In 
addition, studies of the bound states of neutrons in external potential wells meant to 
represent the effects of the neutron-proton interaction may also provide new insights in 
the structure of neutron-rich nuclei. 
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